NYU-TH/98/10/04, HUTP-98/A065, NUB 3187 
BPS Domain Walls in Large N Supersymmetric QCD 



Gia Dvali 1 ' 2 *, Gregory Gabadadze 1 ^ and Zurab Kakushadze 3,4 * 
1 Department of Physics, New York University, 4 Washington Place, New York, NY 10003 
2 International Centre for Theoretical Physics, 34100 Trieste, Italy 
3 Lyman Laboratory of Physics, Harvard University, Cambridge, MA 02138 
^Department of Physics, Northeastern University, Boston, MA 02115 

(January 10, 1999) 



Abstract 

We explicitly construct BPS domain walls interpolating between neighbor- 
ing chirally asymmetric vacua in a model for large N pure supersymmetric 
QCD. The BPS equations for the corresponding Zjv symmetric order param- 
eter effective Lagrangian reduce to those in the A n Landau-Ginsburg model 
assuming that the higher derivative terms in the Kahler potential are sup- 
pressed in the large N limit. These BPS domain walls, which have vanishing 
width in the large N limit, can be viewed as super membranes embedded in 
a (3+l)-dimensional supersymmetric QCD background. The supermembrane 
couples to a three-form supermultiplet whose components we identify with 
the composite fields of supersymmetric QCD. We also discuss certain aspects 
of chromoelectric flux tubes (open strings) ending on these walls which appear 
to support their interpretation as D-branes. 
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I. INTRODUCTION 



Supersymmetric gauge theories provide opportunities to gain analytic control over cer- 
tain complicated non-perturbative phenomena in (3+l)-dimensional strongly coupled gauge 
models The closest supersymmetric cousin of QCD, Af = 1 supersymmetric SU(N) 

gluodynamics (pure supersymmetric QCD), is believed to be a confining gauge theory with a 
mass gap and chiral symmetry breaking H^H. One can therefore expect that a better under- 
standing of pure supersymmetric QCD (SQCD)Q might shed some light on non-perturbative 
phenomena in conventional Quantum Chromodynamics. 

Chiral symmetry breaking in SQCD is due to gluino condensation. The compos- 
ite operator Tr(A a A a ) = AA acquires a non-zero vacuum expectation value (VEV) |Ll§|: 
(AA)fe = N exp(2Trik/N)A 3 , where A is the dynamically generated scale of SQCD. Thus, 
there are N inequivalent vacua related by Z^ discrete symmetry and labeled by the phase 
of the gluino condensate k — 0, . . . , N — 1. In each particular vacuum state the discrete Z^ 
symmetry is broken. As a result, there should exist domain walls separating pairs of distinct 
vacua. 

In H| it was shown that the M = 1 supersymmetry (SUSY) algebra admits a nontrivial 
central extension if BPS saturated domain walls are present in the theory. Assuming that 
such BPS walls exist, their tension can be calculated exactly on the basis of the SUSY 
algebra alone 0. Whether the BPS saturated domain walls are indeed present in a given 
model is a dynamical issue. The strong coupling nature of SQCD, however, makes it difficult 
to study this problem in terms of colored fields. In order to circumvent this difficulty, one 
can use an order parameter effective Lagrangian relevant for the chiral symmetry breaking. 
Thus, in the large A" SQCD one could attempt to solve the corresponding BPS equations in 
terms of colorless composite variables (order parameters). This approach, as we will discuss 
later, can be argued to be reliable in the limit of infinitely large N. This is the route we are 
going to follow in the present work. In particular, we will explicitly construct BPS saturated 
domain walls interpolating between neighboring chirally asymmetric vacua in a model for 
large A" SQCD. The width of these walls is vanishing in the large A" limit. We argue that this 
is precisely the property which gives rise to the wall tension consistent with the prediction 
of the SUSY algebra f|. 

Having found the BPS walls in the large A" limit, we discuss more general questions 
related to physics of extended objects in SQCD. In fact, we address the issue whether the 
lowest-spin states of SQCD can be viewed as excitations of such extended objects. Our 
discussion is motivated by the relation of the large A" gauge models and string theories. 



Indeed, in JT0[ it was shown that in the large A" limit the gauge theory diagrams are organized 



in terms of Riemann surfaces with boundaries and handles, and addition of each extra handle 
on the surface corresponds to suppression by a factor of 1/N 2 . This observation leads to 
the hope that the large A" gauge theory might be described by some kind of string theory. 
Then the large A" expansion of gauge theories would be mapped to the string expansion in 
terms of properly weighted world-sheets of various topologies. 



throughout this paper "SQCD" refers to the J\f = 1 supersymmetric SU(N) gauge theory 
without matter, i.e., pure supersymmetric gluodynamics. 
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The first concrete realization of this idea in the context of weakly coupled gauge theories 
[i.e., when the effective gauge coupling A = Ngy M is fixed at a value A < 1) was given in [11 



in the framework of three-dimensional Chern-Simons gauge theory where the boundaries of 
the string world-sheet are "topological" D-branes^]. More recently, in |]13|L1 string expansion 
was shown to precisely reproduce 't Hooft's large N expansion for certain four dimensional 
(super) co nformal gauge theories. 

In strongly coupled gauge theories, however, the story appears to be much more involved. 
In particular, it is unknown what is the string theory that governs the dynamics of the QCD 
string which is expected to arise as an effective description in the large N limit. This string 
theory is expected to be non-critical flip] , which makes it difficult to study. Nonetheless, in 
|16|1 Witten has made an interesting observation that there might be a connection between 



SQCD strings and domain walls in the large N limit. In particular, based on the assumption 
that the SQCD domain walls are BPS saturated @, Witten argued that in the large N limit 
the domain walls connecting two vacua labeled by k and k! = k + 1 appear to be objects that 
look like D-branes |17| on which the SQCD string should be able to end. Such D-brane-like 
domain walls have also appeared in in the context of large N non-supersymmetric QCD. 

We argue in this work that in the large N limit the BPS domain walls with the vanishing 
width can be viewed as supermembranes embedded in an SQCD background. Moreover, 
the component fields of the three-form supermultiplet which couples to the supermembrane 
can be identified in terms of the composite fields of SQCD. We also point out that some 
of the lowest-spin states of SQCD can be viewed as excitations of a closed string which 
propagates in the vacua separated by the supermembranes. A BPS supermembrane can be 
regarded as a D-brane in the large N limit if there are open strings ending on it. We discuss 



certain aspects of chromoelectric flux tubes ending on the domain walls [19] which appear 
to support the D-brane interpretation of these objects. 

The rest of this paper is organized as follows. In section II we discuss the large N 
behavior of a generic SQCD configuration and point out why a BPS state can exist in 
the large N limit. Then, we briefly discuss the Zn symmetric order parameter effective 



Lagrangian for SQCD (see |20],2lJ and ||22||), and also give a general discussion of domain 
walls in M = 1 supersymmetric theories. In section III we explicitly construct BPS domain 
walls using a model Zjy symmetric order parameter effective superpotential by solving 
the corresponding BPS equations in the large N limit. Note that large iV BPS domain walls 
have been constructed in 



23 



The idea of |23[ is to avoid certain complications arising 
in SQCD by taking an indirect route. In it was argued that in the large N limit 
the problem effectively reduces to the question of finding BPS solitons in the large N A N 
Landau- Ginsburg theory. In this paper we construct the BPS domain walls directly within 
the framework of the Zn symmetric model of large N SQCD. We show that the corresponding 
system of BPS equations indeed reduces to that in the large N An Landau- Ginsburg model 
provided that higher derivative terms in the Kahler potential are suppressed in the large N 



limit. This way we derive the results of 23] which were obtained via an indirect construction. 



2 For recent developments in these directions, see 12]. 

3 This was subsequently generalized to include unoriented world-sheets in [14| 
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In section IV we elucidate some relations to other works on the subject. In section V we 
present our discussions and conclusions. 



II. Zjv SYMMETRIC VACUA AND DOMAIN WALLS 

In this section we review some facts about the large N limit of M = 1 supersymmetric 
QCD. We argue that the SQCD domain wall configurations can saturate the BPS bound if 
the width of the walls vanishes in the large N limit. We will also discuss the symmetric 
order parameter effective superpotential of [^2] which effectively encodes dynamics relevant 
for describing large N domain walls. 



A. Gaugino Condensation and Domain Walls in SQCD 

Consider M = 1 supersymmetric QCD with SU(N) gauge group and no matter. Let A 
be the dynamically generated scale of the model. In this theory the gluino condensate is 
given by g 

(XX) k = Nexp(2mk/N)A 3 . (1) 

Here k — 0, . . . , N — 1, that is, there are N inequivalent vacua corresponding to N different 
phases for the gaugino condensate. Note that the overall factor of N in ([I]) follows from 
the fact that the gaugino condensate (AA) involves a trace over the gauge indices, and the 
resulting VEV should contain the corresponding second Casimir factor. 

According to (H), there are N inequivalent vacua with spontaneously broken discrete 
symmetry in SQCD. Hence, there should exist domain walls separating these different vacua 
H0. On dimensional grounds we expect the tension of such a wall to be of order A 3 . 
It is, however, possible to compute the tension exactly in terms of the gluino condensate 
provided that they are BPS saturated ||. Consider the domain wall separating the vacua 
with (XX) k = Nex.p(27rik/N)A 3 and (XX) v = N exp(2mk' /N)A 3 . The tension of the wall 
(provided that it is BPS saturated) is defined by the central charge Q kk < in the corresponding 
central extension of the M = 1 superalgebra. The central charge is proportional to the 
absolute value of the difference between the values of the superpotential in the two vacua: 

Qkk' oc \W k -W k >\. (2) 

We therefore have the following tension for the BPS domain walls: 

'Tt(k-k')" 



N 2 

1 uui — — 

kk 4vr 2 



sin 



N 



A . (3) 



It is, however, difficult to explicitly construct these domain walls (or even check that they 
are BPS saturated) in SQCD: the order parameter AA is a composite operator, and we are 



A number of interesting properties of domain walls in various models were discussed in [24- 
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dealing with a strong coupling regime. Instead, we will be able to explicitly find BPS domain 
walls in large N SQCD using the symmetric order parameter effective superpotential 



[22 



The presence of BPS saturated domain walls in large iV SQCD has important implications 
|T6|j . Strongly coupled SQCD in the large N limit is believed to be described by a string 
theory with the string coupling X s ~ l/N. Extended solitons in this string theory are 
expected to have tension which goes as 1/A 2 ~ N 2 . Let us, however, consider BPS domain 
walls with k' — k + 1. In the large N limit their tension goes as ~ N ~ 1/A S . In ]16|] it was 
suggested that such domain walls can be viewed as D-branes |T7]] rather than solitons in the 



SQCD string context. Open SQCD strings then can end on these D-branes [16 



B. Domain Walls at Large iV 

Let us consider the domain walls with k' = k + 1 in the large N limit. According to (|3|) 
the tension of such a domain wall (provided that it is BPS saturated) is given by 

T ° = I A 3 - W 

The fact that Tp scales as N in the large N limit might appear a bit puzzling if one thinks of 
such BPS domain walls as SQCD solitons. Indeed, the number of degrees of freedom (gluons 
and gluinos) in SQCD scales as iV 2 in the large A" limit, so naively one expects the energy 
density of a generic SQCD configuration to scale not as N but as iV 2 . If this argument were 
precise, it would imply that the large iV SQCD domain walls with k' = k + 1 could not 
be BPS saturated. However, in the following we will argue that the A" dependence of the 
energy density Tu given by (|j) is indeed correct for the corresponding BPS saturated SQCD 
solitons^. 

In the large N limit the effective Lagrangian is expected to have the following form: 

£ = iV 2 £ (<£,V<l>,V 2 $,...) , (5) 

where $ collectively denotes composite colorless bound states of gluons and/or gluinos. The 
corresponding fields $ are rescaled in such a way that all the leading N dependence in 
(H) is in the overall factor of iV 2 . In particular, the equations of motions for $ do not 
contain any factors of N. Thus, in the large N limit the volume energy density of any 
solitonic configuration e ~ iV 2 A 4 . This implies that the tension of a domain wall T ~ ed is 
proportional to iV 2 A 3 if the "width" of the domain wall d ~ 1/A. However, if d ~ l/N A, 
then we have T ~ NA 3 . This indicates that the domain walls with k' — k + 1 should have 
vanishingly small width d ~ l/N A in the large limit P5[ . In section III we will explicitly 
show that this is indeed the case. Here, however, we will give a simple explanation of why 
d ~ l/N A for such domain walls. 

Since all the fields $ in £ are such that [$[ ~ 1 in the large A^ limit, solitonic solutions 
with d depending on Af might at first seem impossible to accommodate. Thus, naively one 



5 We are grateful to M. Shifman for useful conversations on the issues discussed in this subsection. 
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might expect that Co has no "knowledge" of N. This, however, is not the case provided 
that the number of inequivalent vacua described by £ scales as N []19| , which is precisely 
what happens in SQCD. This way the Lagrangian £ does "know" about N. As a simple 
toy example of such a system consider a single chiral superfield $ with the scalar potential 
V($) = iV 2 |l — & N \ 2 . There are N inequivalent vacua in this theory: = exp(2nik/N), 
k = 0, 1, . . . , N — 1. Thus, even though |$fc| = 1 in V, we have \$k+i — $fc| oc 1/N in the 
large N limit. 

In the following we will show that the model effective Lagrangian for supersymmetric 
gluodynamics [|Z| possesses precisely the properties discussed above. In fact, among other 
things it describes N inequivalent vacua with the broken chiral symmetry. In the large N 
limit it admits BPS solitons with the width proportional to 1/N. These solitons are the 
BPS domain walls separating neighboring chirally asymmetric vacua of the theory. 



C. The Order Parameter Effective Action 

The N chirally asymmetric vacua in SQCD are defined by the VEV of the order parameter 
AA, the gluino bilinear. To construct the domain walls interpolating between different vacua 
we can try to write down an effective Lagrangian for this order parameter. The full effective 
Lagrangian ([5]) would contain infinite number of fields corresponding to higher dimensional 
operators of the theory and higher derivatives as well. However, in the large N limit we 
expect certain simplifications, in fact, the higher dimensional operator contributions should 
be suppressed by extra powers of N. More precisely, we expect this to be the case for domain 
walls interpolating between neighboring vacua; in this case the relative change in the gluino 
condensate (AA) is 0(1/ N), so that we expect the adequate description in the large N limit 
to be in terms of a truncated^ effective Lagrangian containing a finite number of fields. On 
the other hand, if we consider domain walls at finite N the relative change in the gaugino 
condensate (AA) is no longer small, so that truncating the full effective action (§]) to a finite 
number of lowest dimensional operators may no longer be justified. A similar remark applies 
to the case of large N domain walls with \k' — k\ ~ N. In the following, therefore, we will 
focus on the large N limit where we will confine our attention to the domain walls with 
k' = k + 1. 

Thus, we are going to be interested in deducing the effective Lagrangian adequate for 



studying large N domain walls. It is reasonable to assume [20| that the order parameter 



effective Lagrangian should be constructed in terms of the chiral superfield 

S = (Tr(W a W a )) = (Tr(A a A a )) + . . . = (AA) + . . . , 

where S is regarded as a classical superfield, and the matrix elements are defined in the 
presence of an appropriate background source (super)field (for detail see P^| ). Note that 
the gluino bilinear AA is the lowest component of the chiral superfield Tr(W a W a ). (Here 



6 By "truncated" we mean that the corresponding effective action is obtained from the full effective 
action (|5|) by keeping a finite number of fields and neglecting the rest of the fields in the large 
N limit. 
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W a is the usual gauge field strength chiral superfield.) Then, the order parameter effective 



superpotential reproducing all the anomalies of the model is given by [20]: 

W VY = NS 



(6) 



The corresponding scalar potential describes spontaneous chiral symmetry breaking with a 
non-zero gaugino condensate |20| . 



The above effective superpotential, however, is not adequate for studying domain walls 
in SQCD. This can be seen as follows. First, the superpotential (0) does not respect the Z N 



discrete symmetry pl[ . This, in particular, implies that the use of (|6|) cannot be justified 
when describing domain walls interpolating between different chirally asymmetric vacua of 
the model [^1 ] . Thus, the superpotential @ is only adequate for describing a given chirally 



asymmetric vacuum. In fact, as it stands in (g), such a description is only appropriate for 
the vacuum with k = 0. This is related to the strong CP violation. Let us notice that there 
is no CP violation in SQCD since the CP odd 9 angle can always be canceled by chiral 
transformations of massless gluino fields. Let us go back to the effective Lagrangian and see 
how this property is realized there. As we mentioned above, the lowest component of the 
chiral superfield S is the gluino bilinear (AA). The highest component of the superfield S is 
formed by the composite "glueballs": {G^ U G^ U + iG^ v G^ v ). As usual, these latter fields are 
auxiliary components of the chiral superfield S, and, therefore, do not appear in the on-shell 
scalar potential of the modelQ. Thus, the effective scalar potential of the model is a function 
of the gluino bilinear only. Let us now suppose that we would like to deal with the vacuum 
state of the model with a nonzero value of the imaginary phase of the gluino condensate. 
For such a vacuum we have 

Im((AA) fc )=iVA 3 sin^^0. (7) 

The imaginary part of the gluino bilinear is a CP odd quantity. The non-zero vacuum 
value of this quantity signals a strong CP violation in the effective potential. This CP odd 
phase cannot be rotated away as there are no massless fermions in the effective Lagrangian. 
Thus, if the scalar potential is written in terms of the gluino bilinear only, then it cannot 
simultaneously describe the N different vacua and yet preserve the strong CP invariance 
of the model. We must therefore generalize the effective superpotential ([]) by including 
additional fields such that the resulting effective superpotential respects the Z N discrete 
symmetry, and is also consistent with the absence of the strong CP violation. This can be 
done by introducing the "glueball" order parameter fields in the effective Lagrangian. Thus, 



7 These auxiliary fields can appear in the effective Lagrangian if higher super derivatives are present 
in the corresponding Kahler potential. In this case the F-terms become dynamical fields. However, 
including these fields as dynamical degrees of freedom is accompanied by problems with positive 
deflniteness of the corresponding scalar potential (for details see (3(|), and, therefore, it is assumed 
that the Kahler potential does not contain superderivatives. This is expected to be the case in the 
large N limit 
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the observable effects of the CP odd phase in the gluino condensate should be compensated 
by another CP odd quantity, the non-zero vacuum value of the field G liV G ilv . Moreover, this 
should also restore the discrete Z N symmetry. 



The Zjv symmetric effective superpotential has recently been proposed in [ffiql . This 
superpotential involves two chiral superfields 5* and X, where S is the same as before, 
whereas X is a new superfield whose lowest component is related (in a rather complicated 
way) to the gluon bilinears (GG) and (GG) [^,^|. The corresponding scalar potential 
contains the gluino bilinear via the lowest component of the S chiral superfield, and also 
"glueballs" via the lowest component of the X superfield. Indeed, the lowest component 
of X on mass-shell is related to the highest component of S, i.e., to the gluonic operators 
discussed above. Likewise, the F-component of X is related to the lowest component of S, 
i.e., to the gluino bilinear. The explicit relations between these fields are governed by an 
unknown Kahler potential in the on-shell formulation of the model. However, one can think 
of X as a superfield which introduces the necessary order parameters into the description 
and restores the Z^r symmetry of the superpotential @. The effective superpotential for 
the S and X superfields is given by 



W e ff = NS 



(8) 



Here the coefficients c n are 0(1) in the large N limit, and satisfy the following constraints^: 
J2n c n = 0, and Yjn nc n = 1- The additional (compared with (|])) terms in (H) containing the 
X superfield restore the discrete Zn symmetry as well as CP invariance of the effective action 
j22flP|. In particular, the Z^r transformations read: S — > S exp(2ml/N), X — > X — 2ml/ N, 



where I is an integer. Thus, the superpotential (H) describes the N chirally asymmetric 
vacua in SQCD[^. These vacua are given by (S, X) = (Sk,Xk), where Sk = NA 3 exp(— Xk), 
and X k = 2mk/N, k = 0,1, N — 1. 

Here we would like to point out that the effective superpotential (|8]) can be regarded as 
the Veneziano-Yankielowicz superpotential (|5p where the parameter A 3 has been promoted 
to the X dependent chiral superfield A 3 (X) p2| : 

A 3 (*)=A 3 exp(-* + l]r Cn e^ . (9) 



8 The coefficients c n are not uniquely determined within the effective Lagrangian approach, but 
should in principle be fixed via the instanton calculus. If so, it is reasonable to expect that the 
sum over the integer n in (|8|) should contain either only non-negative or only non-positive values 
(albeit this cannot be seen solely on the grounds of the effective Lagrangian arguments). Without 
loss of generality we can then assume that n < in (||). 

9 The corresponding effective action is Z# symmetric as long as the Kahler potential depends on 
X and X* via Z^v symmetric combinations, such as X + X* , e nNX , e nNX * , and so on. 

10 Note that in the corresponding on-shell scalar potential ambiguities related to the different 
branches of the logarithmic function are canceled as well. 



S 



Then the dynamics of the X field defines the phase of A 3 (X) (or, equivalently, the phase of 
the gaugino condensate) in accordance with (2), and the superpotential (§) can be rewritten 
in the following useful form: 



eff 



NS 



In 



+ 1 



The following notation will prove convenient in the subsequent discussions. Let 

X = Ne~ x . 
can now be written as 
A 3 



The effective superpotential 

Weff = NS 

Here f{X) is given by 

/(X) = Xexp 



In 



+ 1 



(10) 



(ii) 



(12) 



1>T (*Y N 
N^ n C ' n \N) 



n>0 



and has the following properties: 

f'(X k ) = , f(X k ) =X k = Nexp(2mk/N) 
Note that f(X) can be rewritten as a polynomial: 



hx)=ny: 



d r , 



n>0 



nN + 



nN+l 



where the sum is over non-negative integers, and the coefficients d r , 
and, result, satisfy the following constraints: 



X d * 







n>0 



X 

n>0 



d r . 



nN + l 



(13) 
(14) 
(15) 

1 are related to c n 's 
(16) 



In particular, d = exp(c /N) = 1 + 0(1/N). 

It is not difficult to see that the superpotential (|T2|) has Z jv discrete symmetry (with the 
transformations S — > S exp(2nil / N) , X — > X exp(2ml/N)), and the corresponding vacua 
are given by 



S — Sh 



X h A 3 



X = X k , k — 0,1, 



N 



(17) 



Thus, the superpotential (|T2|) can be used to describe large N BPS domain walls inter- 
polating between two different chirally asymmetric vacua (S,X) = (Sk,Xk) and (S, X) = 
(Sk>,Xkt) with k' — k + 1. Since the effective superpotential (0) is a truncation of the full 
effective Lagrangian (which unlike ( |T2"D contains an infinite number of fields), it is important 
to understand the applicability limits of this approach. Let the center of a domain wall so- 
lution (which depends only on the z coordinate and is independent of the x, y coordinates) 
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be located at a point z = z . Then, the effective Lagrangian approach should be applicable 
to study the domain wall at distances 

\z - zq\ > d oc , (18) 

where d stands for the "width" of the wall which is vanishing in the large N limit. This is 
the key simplification which allows us to explicitly construct BPS domain walls in the large 
N SQCD. 



D. BPS Domain Walls in N = 1 Theories 

Before we turn to finding BPS domain walls in the large N SQCD, in this subsection 
we would like to review some useful details concerning generic properties of domain walls in 
M = 1 supersymmetric theories. 

Consider an M = 1 supersymmetric theory with one chiral superfield X. Let us assume 
that the corresponding Kahler potential does not contain superderivatives. This is a reason- 
able assumption for an effective theory in the large N limit since higher derivative terms are 
expected to be suppressed by powers oil/N. Let W(X) be the superpotential in this theory 
such that the F-flatness condition Wx = has a discrete set of non-degenerate solutions 
X = X a ,a = l,...,N. 

In such a theory we expect presence of domain walls separating inequivalent vacua X = 
X a and X = Xj,, a ^ b. Let z be the spatial coordinate transverse to the wall. Then 
asymptotically we have X(z) — > X a as z — > — oo, and X(z) — > X b as z — > +oo. Such domain 
walls may or may not be BPS saturated. 

The problem of finding BPS domain walls in four dimensional M = 1 supersymmetric 
theories is really a two dimensional problem as the coordinates x, y along the wall play no 
role in the discussion. Thus, our problem effectively reduces to that of finding BPS solitons 
in two dimensional massive Af = 2 quantum field theories [38|. The corresponding BPS 



equation for the lowest component of the X superfield takes the following form: 

d z X*(z) = e*F* x , 

where F x denotes the F-component of the X superfield, and 7 is some constant phase (whose 
precise value will be given in a moment). If the Kahler potential contains no superderivatives, 
then F x is related to the derivative of the superpotential via the following equation of motion: 
gF x = Wx, where g = g(X,X*) is the Kahler metric defined through the Kahler potential 
K as g = ICxx*- Combining the BPS equation and the equation of motion for F x one 
finds: gd z X* = exp(i-f)Wx ■ Let us now turn to the calculation of the surface energy density 
of a domain wall which satisfies this relation. For a general supersymmetric theory with 
one chiral superfield and with no superderivatives in the Kahler potential the energy of a 
time-independent configuration is given by 

1 r+00 , . 
E ab = - dz(g\d z X\ 2 + V) , (19) 

Z J— 00 ^ ' 

where V = g~ 1 \Wx\ 2 is the scalar potential of the theory. Note that E ab can be rewritten 
as follows: 
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E ab = - dzg- 1 \gd z X* - exp(* 7 )>Vx| 2 + Re (exp(* 7 )(W b - W a )) , (20) 

Z J —CO 

where exp(ry) is an arbitrary constant phase. Note that E ab is independent of 7. For the 
choice 7 = 7 b a we obtain the bound E ab > \W b — W a \ which is the BPS bound. Here 

expH7fea) = W^m • (21) 

The BPS solutions (for which E ab = \Wb — W |) are those that satisfy the following equation: 

gd z X* = exp(z 76a )Wx (22) 

subject to the boundary conditions X(z) — > X a b as z — > =Foo. Going back to the domain 
walls in four dimensions, we have the exact same BPS equation (|22|) , and the tension of the 
domain wall is given by 

T ab = -^\W a -W b \. (23) 

Typically, one does not know the exact form of the Kahler metric g. However, for the 
case of a single superfield X we are considering here one can still make certain statements 
about the corresponding soliton solutions. Thus, let X be a function of a new coordinate z' 
such that it satisfies the following equation 

d z> X* = exp(i lba )W x (24) 

with the boundary conditions X(z') — > X a b as z' — > =1=00. Suppose we are able to find the 
corresponding solution for X(z'). Next, consider the following change of variables: 

d z! z(z')=g(X(z'),X*(z')) , (25) 

where X(z') is the corresponding solution. Note that if we express X as a function of z, 
then X(z) will satisfy the original BPS equation (j22D with the corresponding boundary 
conditions. The change of variables (p5|) is simply a diffeomorphism which is one-to-one as 
long as the Kahler metric g is non-singular. Throughout this paper we will assume that 
the corresponding Kahler metric g is indeed non-singular. Moreover, we will always work 
in the coordinate system parametrized by z', but we will drop the prime for the sake of 
simplicity of the corresponding expressions. (This is effectively equivalent to the case where 
g = 1.) Note that as long as the Kahler metric g is non-singular, the solution of ([24]) (with 
the appropriate boundary conditions) implies that the corresponding solution of ( p2|) exists. 
Moreover, this solution is BPS saturated. 

In the next section we will see that the problem of finding BPS domain walls in the large 
iV SQCD is reduced to that of finding BPS solitons in the so called An Landau- Ginzburg 
model. The solitons in Landau- Ginsburg theories have been studied in detail |38j-440|. Many 
of these models are integrable. The general conditions for existence of BPS solutions were 
first formulated in [|(J, where exhaustive studies of related issues were performed. The 
superpotential in the A n Landau-Ginsburg model (in our notations) is given by 

vN+l 

w K A3 (* - ftt) • (26) 

There are N distinct vacua in this model with X = X^ = exp(2irik/N), k = 0, . . . , iV — 1. 
The solitons interpolating between the vacua X^ and Xy are BPS saturated [p9f|. 
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III. LARGE N DOMAIN WALLS 



In this section we explicitly construct BPS domain walls in the large iV pure super- 
symmetric QCD using the model effective superpotential fll2|) . We will be looking for BPS 
solutions interpolating between neighboring vacua with phases k and k' — k + 1. Before we 
turn to the technical details of the corresponding equations, we would like to make a few 
observations which will lead to important simplifications, and will elucidate some features 
of the system we are dealing with. 

To begin with, let us note that the following relation holds for the vacuum state labeled 
by the phase k: S\k = A 3 f(X)\j,. This relation is nothing but the definition of the gluino 
condensate with the corresponding phase set by the vacuum value of the X superfield. It is 
useful to introduce the chiral superfield NS/A 3 f(X). In all the vacua this superfield takes 
the same value equal N. Let us now concentrate on interpolation between a pair of nearest- 
neighboring vacua. In this case the relative change in the gluino bilinear is of order 1/N. 
Hence, the chiral superfield NS/ A 3 f(X) can only deviate from its vacuum value (which 
equals N) by a quantity of order 1. Thus, one can introduce the following parametrization: 



A3/(X) V N 

with E being a new chiral superfield. Thus, the S superfield is defined via the X and E 
fields as far as the nearest-neighboring vacuum interpolation is concerned. Note that both 
real and imaginary parts of E are at most of order 1 in the large N limit. Substituting the 
expression (p7|) into the superpotential (|3)> one finds 



W eff = NA 3 f{X) 



E 2 / 1 



This is the superpotential describing the variation of the order parameter during the tran- 
sition between nearest-neighboring vacua. The superfield E enters this superpotential in a 
subleading order in 1/N. Thus, it is the superfield X which is left in the leading order, 
and which should describe domain walls between the adjacent vacua in the large iV limit. 
Neglecting higher order terms, the superpotential is therefore given by: 

W eff ~NA 3 f(X) . (28) 

Before we go any further, notice that this description is valid only for a pair of adjacent 
vacua, and a priori the approximations made in (^q ) cannot be justified for a pair of vacua 
for which \k — k'\ ~ N. 

Let us now turn to the Kahler potential describing transitions between a pair of nearest- 
neighboring vacua. Generically, the Kahler potential is a function of the entire tower of 
SQCD composite fields and their superderivatives. It is reasonable to expect, however, that 
in the large N limit the superderivative terms are suppressed by extra powers of 1/N. Thus, 
we assume the following form for the Kahler potential K = Kfi?i,i$f*,X,X*). Here, 
denote all the fields which are generically present in the theory, but which do not appear 
in the large N effective superpotential (p8|) describing the nearest-neighboring transitions. 
Notice that this set of fields includes the superfield E as well (in accordance with (|27|) one 
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of can be denned as = N — £). Here we are normalizing all the fields so that in 
the large N limit they scale as ~ iV (just as is the case for S and X). This, in particular, 
implies that the corresponding components of the Kahler metric are (at most) of order 1 
in the large N limit. Using the Kahler potential and the superpotential ( ^8|) one can write 
down both the equations of motion and the BPS equations. Let us start with the equations 
of motion first. For simplicity of presentation we write down only some part of the system 
of equations of motion relevant for our purposes here. We must make sure, however, that 
eventually all the equations of motion are satisfied. Since there are no superderivatives in 
K, one can write the following equations of motion in the leading order in 1/N: 



where the sign "~" denotes equality in the leading order of the large N expansion. We 
are going to argue below that the fields are irrelevant for the description of the domain 
walls interpolating between adjacent vacua. In other words, we will argue that the domain 
walls between the neighboring vacua are described by the superfield X alone. To see how 
this comes about let us recall that the discrete Z^r symmetry is spontaneously broken in 
the theory. Thus, there are domain walls separating adjacent vacua. The question we are 
dealing with is whether these walls are BPS saturated or not. In other words, there should 
always exist a wall configuration which would satisfy the equations of motion. However, 
the very same configuration might or might not satisfy the corresponding BPS equations. 
Let us start with the equations of motion defined in Given the boundary conditions 

d^(±oo) = and dX(±oo) = 0, and also taking into account non-singularity of the Kahler 
metric, one can rewrite the first equation in (^) as follows: 



While making a transition from one vacuum state to the neighboring one along the solution, 
the relative change in the fields and X should be at most of order 1/N. We therefore will 
be looking for solutions of the form = \l/jo + ^i, where \l/jo (which, if non- vanishing, scale 
as ~ N) denote some constants and set the coordinate dependence of the fields in the 
nearest-neighbor transitions. For instance, as we have already mentioned, the E field enters 
the equations via ^ = N(l — Ti/N). Likewise, the large N behavior of the X superfield is 
set by the expression X — N(l — C/N), where both the real and imaginary parts of ( are 
at most of order 1. Given these relations the coefficients of the Kahler metric, which also 
should be varying by an amount proportional to 1/N, can be written in the following form: 



where CV,- and Cj are some constants of (at most) order 1. Substituting these relations into 
([30]) and taking into account the boundary conditions for the X superfield as well as those 



+ d(gy iX *dX*) ~ , 

3 J 

gxx*F* x + g xn F*. ~ W x , 



(29) 



g* i 9.d** j +gv iX .dX*~0. 



(30) 




(31) 
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for the fields, \l/-(±oo) = 0, one finds that fl30D admits no solutions but the trivial one: 
^ = provided that gm^* | solution oc 0((/N), that is, Ci ~ 0. Note that since domain walls 
must exists (regardless of whether they are BPS saturated or not), the last condition on the 
Kahler potential must be satisfied. On the other hand, as we will show in a moment, once 
this condition is satisfied, the corresponding walls do saturate the BPS bound. Next, to solve 



for C, or equivalently X, one should use the rest of the equations in (|29|) . Let us turn to the 
second equation in (|2D|). As we have shown above the off-diagonal components of the metric 
g^tX* should be of the subleading order on the solution. Thus, in the leading order of the 
large N expansion the second equation in turns into the relation g^^+F^ ,\ so i u tion — 0. 
Consider the trivial solution F^. = 0. This is precisely the solution we would like to consider 
for the purpose of finding BPS domain walls. Indeed, this solution implies that the set of 
BPS equations for ^j, namely, d^i oc F^ i is automatically satisfied (recall that d^i = as 
a result of the first equation in (|29D). 

Finally, let us discuss the third equation in (|2"D|). Given the results of the first equation 
one finds in the leading order: 

g XX *F x ~ W x • (32) 

Let us now turn to the corresponding BPS equation. In accordance with our discussions 
presented above the only nontrivial BPS equation is that for X, d z X*(z) = exp(i7 fc+ljfc )F£, 
where the phase jk+i,k was defined in the previous section. As the next step, one can use 
the equation of motion (^) to derive: 

9xx*d z X* ~ exp(ry fc+lifc )W x • (33) 

Below, we will assume that the corresponding Kahler metric gxx* is non-singular. Then, as 
we reviewed in the previous section, we can absorb the Kahler metric into a redefinition of 
z. (More precisely, we will absorb gxx* I X 2 which is dimensionless.) In fact, in the large N 
limit this procedure is equivalent to the absorption of a finite constant in the redefinition 
of the variable z, or, equivalently, to the redefinition of the width of the domain wall by a 
finite number. The corresponding BPS equation reads: 

d z X* ~ exp(z 7fe+life )r/'(X) , (34) 

where we have introduced the inverse width 

r = dT 1 = NA . (35) 

The boundary conditions for X(z) are given by: X(— oo) = X k , A(+oo) = A fc+1 . 

Thus, the system of BPS equations and equations of motion for this particular case 
reduce formally to the system described by the superpotential 

W ~ iVA 3 /(A) (36) 

amended by the S field defining condition 

S = A 3 f{X) (1 + 0(l/iV 2 )) . 

In addition, the values of the all other fields ^[ (including S) on the solution vanish in the 
leading order of the 1/N expansion. These relations are supposed to be satisfied on the 
particular BPS solution we are dealing with. In the following subsections we will discuss the 
solution to the BPS equation (|34]) and present the explicit form of the corresponding BPS 
domain wall. 
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A. A Simple Example 



Before solving the BPS equation (|34D for the most general form of f{X), in this subsec- 
tion, for illustrative purposes, we will consider a simple example where we take the values 
of the numerical coefficients d n as follows: d$ = —d\ = 1 + l/N, and all the other d n (n > 1) 
coefficients are 0(1/N). Then in the leading order in l/N one finds: 

N /X\ N+1 

, W=! *(*-_(-) )• (37) 

Note that the effective superpotential (^) in this case reduces to that of the model 
discussed in the previous section. 

Next, let us solve the BPS equation (|34D in the large N limit for f(X) given by (|37|). (In 
this subsection we will closely follow the corresponding discussion in ||23|| .) Let X = X^Y. 
Then we can rewrite (0) in terms of Y: 

Nd z Y* ~ -zexp(-7ri/iV)r [l - Y N ] . (38) 

The boundary conditions on Y read: 

Y{z -> -oo) = 1 , Y(z -> +oo) = exp(27ri/iV) . (39) 

Thus, in the large N limit arg(F) changes by a small amount, namely, from to 2tc/N. It 
is convenient to parametrize Y as follows: 

Y = (1 - p/N)exp(i((f) + 7r)/N) . (40) 

Here p and 4> are real. The boundary conditions then read: 

p(z -> ±oo) = , (f)(z -> ±oo) = ±7T . (41) 

Then (|38|) becomes a system of the following first order differential equations: 

d z <j) = T [1 + exp(-p) cos(0)] , (42) 
d z p = — rexp(— p) sin(0) . (43) 

Here we are taking the large iV limit and only keeping the leading terms. In particular, we 
have taken into account that (1 — p/N) N — > exp(— p) as N — ► oo. 

The above system of differential equations can be integrated. Note that these equations 
do not explicitly contain z. This implies that if <p(z) and p(z) give a solution with the 
appropriate boundary conditions, so will §{z — Zq) and p(z — zo) for any constant Zq. (This 
is simply the statement that the system possesses translational invariance in the z direction.) 
The general solution is given by 

cos(0) = (p — 1) exp(p) , (44) 

[ P{Z) di [exp(-20 - (1 - i) 2 Y l = -Y\z - z \ , (45) 

Jp(z ) 

where p(z ) = po(~ 1-278) is the solution of the following equation: 
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(po - 1) exp(po) = 1 



(46) 



Note that z is the center of the domain wall. These solutions have characteristic for a 
domain wall behavior. The amplitude p changes from zero at — oo to zero at +oo with a 
bell-shaped extremum at z = zq. Likewise, the phase <fi changes monotonically from —it at 
— oo to +7r at +oo going through the origin at z = Zq. The width of this wall d ~ 1/NA 
vanishes in the large N limit. 



B. The General Case 

In this subsection we construct the BPS domain walls for the most general form of 
f(X). In fact, as we will see in a moment, the corresponding BPS equations are actually 
reduced to those discussed in the previous subsection. That is, for the most general effective 
superpotential W in ([36]) the problem of the large N domain walls can be reduced to that 
of finding BPS solitons in the large N Landau- Ginsburg theory. 

To begin with, let us simplify the form of the effective superpotential as follows. Let us 
make the following change of variables: 

m= m=x-jb{ii)™- (47) 

This transformation preserves the Z^ symmetry: in terms of the new variable X we still 
have iV non-degenerate vacua X = X^ = N exp(2iTik/N), k = 0, 1, . . . , N — 1. Also, the 
above non-linear change of variables is non-singular as long as the coefficients d n satisfy 
constraints (|i~6|). This implies that the resulting Kahler metric g(X,X*) is non-singular as 
well (provided that the original Kahler metric g(X,X*) was non-singular). Therefore, as 
we reviewed in the previous section, we can absorb the Kahler metric into a redefinition of 
z. (More precisely, we will absorb g/A 2 which is dimensionless.) The corresponding BPS 
equation reads: 

d z X* = exp(i 7fc+life )rf (X) , (48) 

where the phase jk+i,k is the same as in the previous subsection. The boundary conditions 
for X(z) are given by: X(— oo) = X k , X(+oo) = X fe+1 . Note that the BPS equation 



(|48|) is the same as that derived in [23] using a more indirect construction. In particular, 
it coincides with the BPS equation we just solved in the previous subsection for a simple 
example corresponding to the A N superpotential. Thus, the problem of BPS domain walls 
in the large N SQCD is indeed reduced to that of finding BPS solitons in the large N 
An Landau- Ginsburg theory. In fact, this is a generic property of Zn symmetric Landau- 
Ginsburg theories - they are all related via non-linear change of variables which is non- 
singular. This implies that if the Kahler metric is non-singular for one choice of variables, 
it is non-singular in the transformed variables as well. Since the induced Kahler metric can 
always be absorbed into the redefinition of z, this allows us to solve the BPS equations 
exactly (in the large N limit) in the suitable coordinate "frame" . 

For completeness we note that the shape of the gaugino condensate is given by 
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S = A 3 f(X) = A 3 /(X) ~ A 3 X 



(49) 



Thus, in the large N limit the gaugino condensate is controlled by the VEV of the chiral 
superfield X. This important simplification is the key observation of |p3 |. In particular, 



as we explain in section IV, the additional singlet superfield used in p3[ can be related to 



X. This point is central to the fact that the indirect construction of gives BPS domain 
walls in the large N SQCD, and not in some other theory. 
The tension of the above domain walls is given by: 



T; 



D 



N 



A 



(50) 



The width of these domain walls is ~ r~ x which goes as ~ 1/N. Thus, these domain walls 
are infinitely thin in the large N limit. Moreover, they are BPS saturated. Note that there 
are other domain walls in large N SQCD p3|| , namely, those with k' ^ k ± 1. We refer the 
reader to |23| for details. 



C. Domain Walls and the Three- form Supermultiplet 

The large N SQCD domain walls found in the previous sections can be viewed as mem- 
branes which break half of SUSY generators and are embedded in a (3 + l)-dimensional 
SQCD background. If so, a three-form supermultiplet |?]J must couple to such a membrane 



to have world-volume K-symmetry j42]0. If the domain walls are indeed 2-branes (mem- 
branes), we must be able to identify the corresponding three-form supermultiplet with the 
appropriate (3 + l)-dimensional "bulk" physical fields of SQCD. In fact, in the following 
we will argue that after gauge fixing the three-form superfield reduces to component fields, 
which in a given chirally asymmetric vacuum describe the corresponding gluon and/or gluino 
composites of SQCD. 

Before we turn to this identification, let us make the following remark. Consider the 
domain wall separating two neighboring chirally asymmetric vacua. In each of these vacua 
there are SQCD fields with a well defined particle interpretation in terms of colorless com- 
posite bound states of gluons and/or gluinos. We would like to show that the SQCD fields 
in each of the two vacua consistently couple to the membrane world-volume. Note that the 
corresponding couplings for these two vacua are different. A consistent "sewing" of these 
couplings is guaranteed by the fact that the domain wall is a solution (interpolating between 
these two vacua) of the order parameter effective action (|8]). In the following we will there- 
fore identify the components of the three-form supermultiplet with SQCD colorless fields in 
a given chirally asymmetric vacuum. 

Let us start with the supersymmetric Green-Schwarz action for a 2-brane (membrane) 
embedded in the (3 + l)-dimensional target space. The membrane action S m includes the 
Wess-Zumino term corresponding to the three-form coupling to the world-volume of the 
membrane: 



11 As usual, K-symmetry eliminates half of the world-volume fermionic degrees of freedom, which 
is necessary for having a supersymmetric formulation of the membrane action. 
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S m = T D J d 3 a (--^gg^ElEi Vij + -V=S + -e abc E*E*ECB ABC ^ , (51) 

where a a (a = 0, 1, 2) denote the membrane world- volume coordinates, and Q a b is the induced 
metric on the world-volume. Here the pull-back E^ = 8 a Z M {a)Ey [ i Z M = (x^, 9 a ) are the 
target superspace coordinates, and E^ denotes the supervielbein, where A = (i, a) are the 
tangent superspace indices. The constraints pfl||i3|| imposed by K-symmetry on the three- 
form superfield Babc can be solved, and the remaining independent components of Babc 
can be combined into the following real tensor superfield: 

U = B + i6x-i9x+ —0 2 A* + —6 2 A + —9a l "Be uvXp C ljXp + 

16 16 48 M p 

\<?0 ^ + 0*8^ + \d 2 - + -/0 2 Qe - d 2 B) . (52) 

We would like to point out that the supermembrane action by itself gives rise to the real 
tensor superfield fl52|) only defined up to a shift by a linear supermultiplet L satisfying 
D 2 L = D 2 L = |0],[E^]. If the membrane is considered in empty space, then this shift, 
U —>■ U — L, can always be used to remove half of the components in (0). However, in 
our case the membrane is embedded in a space where composite colorless SQCD fields live. 
Some of these fields, as we are going to see shortly, couple to the membrane. Thus, the total 
classical action of the system is the sum of the membrane action S m and the effective action 
for the bulk fields 

S = S m + Sbuik ■ (53) 
The action for the bulk fields S^uik is a rather complicated object and generically consists 



of an infinite number of fields. However, it was shown in 137] that lowest-spin SQCD exci 



tations (spin-zero glueballs, gluino-gluino mesons and their fermionic superpartners) can be 



described by the components of the three-form supermultiplet U. It was also argued in |37 
that the shift U —* U — L is not a symmetry of the bulk effective Lagrangian S^uik- Thus, all 
the components in ( |52f) should be retained as physical ones. If so, one can identify the fields 
in (|52| ) with the bulk composite massive fields of SQCD j37j. Thus, the three-form field 
C^x is just the "magnetic" dual of the Chern-Simons current of SQCD, and is related to 
the operator GG via e^ v x p d^C vXp oc GG. (Note that C pv \ is a massive field in this approach 
with one physical pseudoscalar degree of freedom [0.) The two Weyl fermion states x an d 
^ in (p2]) correspond to (l,s) = (0,1/2) and (l,s) = (1,1/2) gluino-gluon bound states, 
respectively (where I and s denote the orbital momentum and spin respectively). The A 
field denotes the AA state, whereas the B field is related to the S field (on-shell) which 
corresponds to the scalar glueball given by the composite operator. 

To summarize, there are four real massive scalar degrees of freedom corresponding to the 
scalar and pseudoscalar gluino-gluino composites, and the scalar and pseudoscalar glueballs. 
In addition, there are four real (two Weyl) fermionic degrees of freedom corresponding 
to the gluino-gluon composites with different orbital momenta. These states are present 
in the superfield fl52]) which necessarily appears in the supersymmetric formulation of the 
membrane action. 

One might wonder now what is the relation between the bulk effective action <S& u z/- and 
the order parameter effective action for two chiral superfields S and X used previously. As we 
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mentioned above, S^uik depends in general on an infinite number of superfields corresponding 
to different spin-orbital excitations of SQCD. If one restricts consideration to the lowest-spin 
states only, then those states precisely fit in the supermultiplet fl52|) pTfl, so that Sbuik — 



<Sbuik(U, ...)• On the other hand, it has been known for some time |4^ , ^5| that the field 
content of the three-form supermultiplet U can be rearranged into two chiral superfields. 
This property was used in 0,^] to rewrite the three-form supermultiplet U in terms of two 
chiral superfields which in our case, in a given chirally asymmetric vacuum, correspond to 
the chiral superfield S and some function of the chiral superfield X ||22|| . 

Before we turn to the next section we would like to comment on SQCD strings and some 
related phenomena. There are a number of indications, as we discussed in the introduction, 
that in the large N limit QCD (or SQCD) can be viewed as some non-critical string theory. 
Thus, one might be able to describe all the SQCD states in the large N limit as string 
excitations, without actually referring to the fundamental theory (SQCD). The question we 
address here is whether some of the SQCD states discussed above could indeed be identified 
with excitations of strings present in the "bulk" . In the string picture glueballs are identified 
with excitations of a closed non-critical string: the scalar glueball in the SQCD string theory 



would correspond to the string "dilaton" jR| . The pseudoscalar glueball is the analog of the 
"axion" and, by supersymmetry, a gluino-gluon bound state corresponds to the "dilatino" . 
These states, as we discussed above, along with the gluino-gluino bound states, are combined 
in one massive tensor supermultiplet, similarly to what can happen in superstring theories 



at the non-perturbative level (see discussions in the second reference in |[45|| ). Note that 
the closed string itself, whose excitations we identify with the glueballs, can be viewed as 
a limit of a closed membrane with the topology of a torus. This is consistent with the 



results of [gg] where a part of the low-energy glueball spectrum of QCD was calculated in 
a closed bosonic membrane model and a good agreement with the lattice QCD predictions 
was foundQ. These observations might provide some hints toward understanding the string 
theory of strong interactions. 

Finally, we would like to comment on the possibility of SQCD open strings ending on 
the domain walls. If the domain walls we found are to be interpreted as D-branes, there 



should be open strings ending on the brane [16]. Thus, in the D-brane picture, one should in 



principle be able to identify the U(l) gauge field in the membrane world- volume action [T7j. 
One might wonder how this U(l) gauge field arises in the field theory context if the open 
strings we discuss are viewed as SQCD chromoelectric flux tubes. An example of a wall on 
which the flux tube can end was constructed in |J. In that case, the corresponding theory 
is in the confining phase outside of the wall. Thus, all chromoelectric fields are squeezed 
into flux tubes outside of the wall. It is energetically favorable for the tubes to end on the 
wall and to spread the flux into the wall interior if the gauge theory inside the wall is in the 
Higgs phase |§ . Thus, the original gauge group of the confining theory, being Higgsed inside 
the defect, could support flux tubes ending on the wall. The Higgs phase inside the defect 
can support the residual unbroken U(l) gauge group in the wall. In the large N limit the 
wall becomes infinitely thin. Thus, the system of the flux tube and the wall can be regarded 



12 In p6fl , because of computational difficulties, the spectrum of glueballs was actually calculated 
for a bosonic membrane with the topology of a sphere. 
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in that limit as an SQCD string ending on a D-brane. The world- volume U(l) field in that 
case can be interpreted as a collective solitonic excitation in the soliton picture for the wall, 
or, alternatively, as a mode of an open string ending on the D-brane, if the D-brane picture 
is adopted. Recently, some features of the system of N such domain walls with flux tubes 
were discussed in |T7| . It would be interesting to study these issues further. 



IV. RELATIONS TO PREVIOUS STUDIES 

Surprisingly enough, the solution of BPS equations as well as superpotential fl37|) coincide 



with the ones derived in [23] within a completely different approach. As we will argue below, 



this coincidence is not an accident. In [23[ BPS domain wall solutions were found in the 
theory obtained by integrating out heavy states in SQCD with Nf heavy quark flavors. The 
method of gaining an analytic control over BPS equations in pure SQCD by adding flavors 
and then integrating them out was used earlier in a number of papers, and, in that respect, 



was not new in p3| . However, there is a crucial point which makes conclusions of [23| so 
different from other approaches. Let us explore this difference in somewhat more detail. As 
is well known, the supersymmetric QCD with Nf < N flavors possesses an anomaly free 
R-symmetry under which flavors carry the charge Nf N N ■ Adding a constant mass term for 
"mesons" M in the supersymmetric Lagrangian 

mTr(M) (54) 

explicitly breaks this symmetry. However, in the massive theory there still is an anomaly 
free discrete Z^v symmetry. This group acts on the gaugino bilinear, and, in fact, coincides 
in that respect with the anomaly free symmetry of pure SQCD. Moreover, spontaneous 
breaking of this discrete symmetry by meson VEVs should produce domain walls which 
after taking the limit m — > oo are expected to transform continuously into the domain walls 
of pure SQCD. This was the method used in f31~f . It was shown in |31]] that within this 
approach the BPS walls do exist for sufficiently small values of m, however, above some 
critical the walls cease to exist. Thus, one would seemingly not have BPS domain walls 
in pure SQCD which is expected to be recovered in the limit m — > oo. Although per se there 
is nothing wrong with solving equations for fixed m and then extrapolating them into the 
large m region, it seems to us that the domain wall solutions obtained this way cannot be 
related to the walls in pure SQCD discussed in M and explicitly constructed in the present 
work. The reason for this is as follows. The heavy mesons which are being integrated out 
in this approach transform under the discrete symmetry which also acts on the degrees of 
freedom of the corresponding low-energy theory. If so, a straightforward decoupling of these 
fields cannot be justified for the problem of finding BPS domain walls of the corresponding 
low-energy theory. Indeed, it was shown by Kogan, Kovner and Shifman in |L9| that when 
the heavy fields are charged under the discrete symmetry (which is subsequently broken 
spontaneously) the classical solutions that interpolate between different vacua cannot be 
adequately described by simply integrating heavy fields out and solving equations of the 
corresponding effective low-energy theory. This inadequacy arises due to "cusp singularities" 
which emerge in the low-energy theory once the fields charged under the discrete symmetry 
are integrated out [I9|. Now, in SQCD with flavors and the regular mass term, mesons 
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necessarily transform under the remaining discrete symmetry group, and, thus, the "cusp 
crossing" will always occur once the mesons are integrated out. In this case, albeit one can 
consistently solve the BPS equations and find the wall solutions for a finite mass m ||31|| , 



this method would not allow one to account for the BPS walls of pure supersymmetric QCD 
which is reached in the m — > oo limit. In order to avoid this difficulty, the mesons which are 
being integrated out should not transform under the discrete symmetry. This is possible to 
accomplish if the mass term is not a constant, but rather is promoted into a dynamical field. 
This was the key observation of p3|j . Consider SQCD with Nf = N flavors. This theory has 
an anomaly free £/(l)i?-symmetry under which mesons are neutral. To avoid the presence 
of "cusp singularities", one can postulate that the mass term for "mesons" arises due to 
the VEV of a new chiral superfield X, and takes the form AXTr(M). In this case there is 
a discrete Zjy symmetry under which the X field and gauginos transform but mesons are 
neutral. This group, in fact, is a subgroup of the anomaly free U(1)r. Now, if the mesons 
acquire large masses due to the VEV of the X superfield, then these mesons can be safely 
integrated out without encountering cusp singularities discussed above. Thus, one should be 
able to describe SQCD domain walls in this case. The necessary condition is that one should 
work in the region where X is large. For this, one can add a small Z^-conserving term in the 
superpotential. The simplest one is X N+1 /(N + 1). In this case, after integrating mesons 
out, one recovers (p7|). Another indication of a deeper connection between X in (|T2| ) and the 
one in [23] is the fact that in the above context X essentially acts as a dilaton (plus axion) 
field. This can be seen in a number of ways. The simplest one is to point out that the VEV 
of X spontaneously breaks an anomalous {7(l)-symmetry (the Peccei-Quinn symmetry p9| ) 
and, thus, its imaginary part is an axion field with a one-loop anomaly-induced coupling 

one loop factor x lm(X)FF . (55) 

The supersymmetric generalization of this coupling gives 

\n(X)W a W a . (56) 

Thus, ln(AT) essentially sets the value of the inverse gauge coupling 1/g 2 . This precisely 
matches (|T2"D and the form of the effective superpotential derived from gaugino condensation 
in the low-energy theory for large X. Indeed, for large X the effective low-energy theory is 
pure Yang-Mills (plus X) with a gaugino condensate. The effective superpotential is simply 
(AA) ~ A 3 where A is a low-energy scale of the theory which from the one-loop matching of 
the gauge couplings at the scale X is simply A 3 ~ VAq CD . This is in agreement with the 
dilaton superpotential induced by the gaugino condensation 

XA' QCB ~ e"^ 2 A 3 CD . (57) 



According to |£8|, the dilaton of QCD is related to Cfr u operator, which gives another way of 
deriving why X should be interpreted as the field responsible for "gluonic" (that is, glueball) 



degrees of freedom in accordance with [22] and our analysis in this paper. 



V. DISCUSSIONS AND CONCLUSIONS 

In this section we briefly summarize the main results of our discussions. Before we turn 
to the conclusions we would like to make the following comments. 



21 



Suppose we consider now domain walls at finite N. Can we use an effective superpotential 
(say, of the type we used here) to study them in this case? It appears that using any type of 
truncated effective superpotential may not be justified. The reason why is that the gluino 
condensate changes substantially inside of the domain walls, so that integrating out heavy 
fields, let us call them Y m , might not be justified. Does this imply that SQCD domain 
walls are not BPS saturated at finite A? Strictly speaking one cannot draw this conclusion 
as the effective description in terms of S and X only may not be valid, and one might 
have to include all of the infinitely many fields X, S, Y m to obtain an adequate description. 
Practically, this means that the problem may not be easily tractable at finite N. 

To summarize, we have constructed the BPS domain walls interpolating between neigh- 
boring vacua with broken chiral symmetry in large N supersymmetric QCD. The tension of 
these walls, saturating the BPS bound, scales as ~ N, and the width of the walls vanishes as 
~ 1/N in the large N limit. These walls can be interpreted as supermembranes embedded in 
a (3+l)-dimensional SQCD background. The components of the three-form supermultiplet 
that couples to the supermembrane were identified in terms of composite fields of SQCD. 
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